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'  A  study  is  being  made  of  the  effect  oh  aerodynamic  drag  of  boundary  layer 
irradiation  by  radioactive  sources.  A  blow-down  wind  tunnel  and  a  skin  friction 
drag  balance  have  .been  designed. and,  constructed.  The  balance  is  of  the  null- 
position  type,  and  is  operated  using  an  automatic  control  system  to  maintain  the 
null  position  of  the  drag  plate.  High  accuracy  and  stability  are  observed  at 
flow  velocities  up  to  200  m/s  and  resolutions  of  up  to  0.1%  of  changes  in  drag 
have  been  achieved.  In  separate  experiments,  a  study  of  the  effect  of  radio¬ 
active  emission  on  gas  viscosity  is  being  made  with  a  specially  designed  torsion 
disc  viscometer.  Measurements  to  date  have  been  at  atmospheric  pressure  with 
both  pieces  of  apparatus.  No  significant  changes  in  drag  have  been  found  except 
at  low  f low  speeds.  However,  future  work  will  be  extended  to  lower  pressures  wher  * 
there  are  reports  of  a  decrease  in  viscosity  with  radiation.  An  objectiv  > 

of  new  research  will  be  to  find  the  optimum  conditions  for  any  drag  reduction.  ^ — 
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Both  forms  of  apparatus  have  the  sensitivity  and  flexibility  to  study  aerodynamic 
drag  and  gas  viscosity  for  a  variety  of  configurations.  ( 
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•i ..  INTRODUCTION'  • 

•  *  '-  “  1  -  11 V"  ,Tl  "  r'  'i 

A  previous  report‘d  described  in  detail  the  development  ob^an— aerodynamic 
'skin  friction  drag  balance  for  drag  reduction  studies:  with  radioactive  irradiation 
of  the  air  flow.  This  report  describes  the  continuation  of  .this  research.  It 
includes,  an  -analysis  of  possible  error  sources,  improvements  in  the  drag, measure¬ 
ment  technique,  measurements  at  lower  flow  speed*  the  measurement  of  skin  friction 
drag  in  both  laminar/ transitional  and  fully  turbulent  flow,  and  an  exhaustive 
set  of  measurements  on  drag  reduction  by  irradiation.  In  addition,  some  research 
using  torsion  disc  viscometry  to  investigate  the  effect  of  irradiation  on  gas  viscos¬ 
ity  has  been  initiated. 

In  view  of  the  detailed,  description  of  the  .development  and  construction  of 
the  drag  balance  given  in  the  previous  report,  only  a  brief  resumd  of  its  salient 
features  will  be  given  in  this  report,  which  will  concentrate  on  the  experimental 
data  obtained,  with  the  balance. 

The  torsion  disc  work  is  at  an  early  stage.  An  account  is  given  of  the 
design  .and  construction  of  a  prototype  torsion  disc  viscometer,  together  with 
an  analysis  of  the  data  obtained  so  far. 

; 

2.  DRAG  BALANCE  EXPERIMENTS  ’  t 

The  drag  balance  is  of  the  floating  element  type  and  is  usually  operated 
in  the  null  position  mode,  the  drag  force  being  balanced  by  coils  and  magnets 
■acting  so  as  to  keep  the  drag  element  in  the  null  position.  The  displacement 
is  measured  very  accurately  using  a  linear  variable  differential  transformer 
(LVDT),  (See  Figs.  1  and  2).  The  balance  can  also  be  operated  in  the  deflection 
mode.  This  is  done  for  very  small  forces  (<  300  mg  force),  the  drag  being 
determined  from  the  LVDT  voltage  corresponding  to  the  deflection  of  the  drag  < 

element  against  its  resilient  piano  wire  suspension. 

The  centre  line  flow  velocity  at  the  vindtunnel  test  section  is  determined 
using  a  miniature  pitot  tube  and  static  orifice  just  downstream  of  the  drag 
element.  The  pressure  is  measured  on  a  precision  dial  manometer  for  pitot  static 


pressures  greater  than  !  mmHg.  For  lower  pressures  a  tilt  paraffin 
manometer  is  used.  The  static  pressure  in  the  balance  housing  is  also 
monitored  with  a  paraffin  manometer  as  a  check  for  leaks  which  would 
affect  the  drag  readings.  The  air  temperature  in  the  housing  is  monitored 
by  thermocouples. 

In  the  null  position  mode,  the  null  position  is  maintained  by  an  automatic, 
control  circuit.  This  is  triggered  by  the  displacement  of  the  plate  as  measured 
by  the  LVDT.  The  coil  current  necessary  to  maintain  the  null  position  is  measured  i 
terms  of  the  voltage  across  a  2.7  55  resistor  in  series  with  the  coils.  This 

voltage  is  fed  to  a  DVM  and  to  an  X-Y  recorder.  The  true  drag  voltage  can 
be  most  easily  determined  from  the  fluctuating  voltage  (due  to  vibration  and 
fluctuating  flow  speed)  by  looking  at  a  trace  on  the  recorder  for  a  period 
of  5  seconds  or  so.  This  method  has  advantages  over  an  integrating  circuit  and 
voltmeter,  as  the  response  of  the  drag  balance  to  a  vibration  of  moderate 
amplitude  is  easily  recognisable,  and  can  be  ignored,  whereas  the  integral 
of  a  damped  simple  harmonic  vibration  would  be  non  zero.  For  deflection 
measurements  at  low  speed  and  high  sensitivity,  an  integrating  circuit  has 
to  be  used  to  eliminate  low  frequency  hum  and  to  protect  against  large 
amplitude  vibrations.  Fig.  3  shows  a  block  diagram  of  the  windtunnel  instrumen¬ 
tation. 

In  order  to  make  drag  measurements  in  fully  turbulent  flow  a  0,91  mm  dia. 

(20  SWG  copper  wire)  boundary  layer  "trip"  can  be  placed  against  the  upper 
windtunnel  surface  at  a  point  110  mm  upstream  of  the  leading  edge  of  the  drag 
plate.  This  gives  fully  turbulent  flow  down  to  the  critical  channel  Reynolds 
number. 

2. 1  Summary  of  Previous  Research 

The  previous  work  showed  that  the  drag  balance  and  control  circuit 
worked  well,  and  that  it  was  possible  to  measure  drag  forces  with  good 
accuracy  and  consistency.  However,  some  problems  were  experienced  with 
the  construction  and  mounting  of  the  radioactive  plates  in  the  drag  balance. 


3. 


Some  tentative  results  we re  obtained  which  showed  a  drag  reduction*  but  the 
accuracy  was  poor  and  it  was  not  possible  to  make  an  extended  series  pi 
readings.  ' 

2.2.  Analysis  of  Error  Sources 

As  mentioned  in  the  previous  .report,  it  is  hot  necessary  to  obtain 

an  absolute  determination  of  the  drag  coefficient.  This  is  because  systematic 

errors  due  to  factors  such  as  plate  misalignment  or  pressure  gradients  only 

resuit  in  a  reduction  of  the  sensitivity  of  the  drag  reduction  measurements. 

For  example,  a  spurious  contribution  of  10%  in  the  drag  force  would  reduce 

u  1%  resolution  in  drag  change  to  1.1%.  However,  the  magnitude  of  some  of 

these  effects  was  investigated  to  ensure  that  they  were  of  a* low  enough  order,. 

Of  primary  importance,  though,  is  the  degree  of  consistency  between 

measurements  with  different  active  plates  mounted  in  the  drag  balance.  Due 

to  the  imperfections  of  the  surfaces  of  the  plates,  and  the  difficulty  of 

getting  a  sufficiently  accurate  alignment  with  the  lower  windtunnel  surface, 

the  drag  measurements  were  not  adequately  reproducible.  It  was  decided  to 

instead  have  the  radioactive  foils  not  in  contact  with  the  flow,  but  to  fix 

a  1/16"  perspex  (PMMA)  shield  in  the  lower  wall  of  the  windtunnel.  This 

gives  a  20%  reduction  in  the  radiation  level  near  the  drag  element  for  the 
147 

Pm  sources,  emitting  0.22  McV  3  particles.  This  is  more  than  compensated 

210 

for  by  the  improved  accuracy  of  measurement.  The  a  source  (  Po)  cannot, 

58 

of  course,  be  used  in  this  configuration;  the  yt  Co)  source  can,  but  these 
experiments  have  not  yet  been  performed.  As  the  sources  are  no  longer  in 
contact  with  the  flow,  the  exhaust  duct  is  not  needed,  permitting  more 
frequent  calibration. 

As  mentioned  at  the  start  of  this  section,  measurements  have  been  made  at 
lower  flow  speeds.  For  low  speed  null  position  measurements,  the  LVDT 
deflection  ''oltage  is  comparable  with  the  trigger  levels  of  the  control  circuit. 


It  was  thought  that  this  might  lead  to  inaccuracies  *  but  in  fact  trigger  ' 

> 

levels  of  two  or  three  times  the  standard  did  not  'signi£icantl.y«rlegrade  the 
performance  of  the  control  circuit.  At  medium  speeds  (Fig.  4)  the  difference 
observed  is  -of  the  order  of  the  consistency  of  the  drag  measurement. 

To  evaluate  the  effect  of  the. pressure  difference  on  the  ends  of  the 

* 

drag  element,  a  special  drag  plate  was  constructed  with  i  mm  dia;  static 
taps  in  the  ends,  at  1 f  and  3  mm  from  the  surface  of  the  plate.,  connecting  with 
small  capillary  paraffin  manometers  suspended  with  the  drag  element.  These 
registered  pressure  differences  of  as  little  as  0.02  mmllg. 

Measurements  made  over  a  range  of  flow  velocities  show  that  the  effect 
of  the  pressure  gradient  is  by  no  means  simple,  in  fact,  there  can  even  be 
a  pressure  force  acting  upstream  (see  Fig.  5;  a  +ve  difference  corresponds 
to  a  downstream  force,  -ve  an  upstream  force).  This  force  is  also  influenced 
by  the  orientation  of  the  plate,  i.e.  the  up  and  downstream  gapsize,  and  more 
strongly,  by  the  alignment  with  the  windtunnel  wall  (see  Fig.  6.  Data  are 
for  the  static  tap  at  1.5  mm). 

Fig.  7  shows  the  typical  %  contribution  of  the  pressure  force  to 
the  total  force.  The  calibration  of  the  manometers  had  to  be  estimated,  and 
may  be  rather  inaccurate,  Several  measurements  taken  at  positions  of  zero 
pressure  difference  show  that  the  basic  character  of  the  drag  curve  is  un¬ 
affected  by  the  pressure  force  However,  as  in  the  recent  series  of  experiments 
the  radioactive  plates  were  mounted  in  the  lower  windtunnel  wall,  and  not  as 
a  drag  element,  it  was  possible  to  redesign  the  drag  element  with  the  ends  out 
at  45°  to  a  sharp  edge.  This,  of  course,  reduces  the  effect  of  the  pressure 
gradient,  although  the  system  is  possibly  more  sensitive  to  misalignment. 

2.3  Experimental  Results 

Before  beginning  measurements  with  the  radioactive  sources,  a 
comprehensive  set  of  drag  measurements  was  Laken  to  check  the  performance 
of  the  system.  Fig.  8  shows  a  typical  calibration  curve  for  the  balance 
in  the  deflection  mode;  fig.  9  a  calibration  curve  for  the  null  position 


shown  in  Figi  .10. 

It  can  be  seen  that  with  the  boundary  layer  trip  in  place  the  flow 
becomes  turbulent  above  a  channel  Reynolds  number  of  about  2500  arid  is. 
essentially  fully  developed  by  about  Re  =  4000.  Without  the  trip  the  flow 
remained  laroiriar  up  to  Re  =  9000,  when  transition  to  fully  turbulent  flow 
beings.  Although  transition  seems  complete  by  Re  =  30,000  it  can  be  seen 
that  the  drag  coefficient  is  still  less  than  before.  This  is  probably  due 
the  turbulence  being 'initiated  at  a  different  position  along  the  windtunnel. 

The  dat$i  are  in  good  agreement  with  ref.  2  and  3. 

Measurements  with  radioactive  sources  were  then  made  in  the  null  position 

mode.  The  measurements  are  presented  in  Tables,!  and  2,  *  denotes  a  measurement 

,  ♦ 

which  has  been  estimated  for  computational  purposes  as  the  original  reading, 
was  either  missing  or  more  than  3  standard  deviations  from  the  mean.  Plates 
N1  to  N5  are  inactive;  A1  and  A2  contain  approximately  100  and  50  mCi  of  *  7 Pro 
respectively. 

The  measurements  have  been  corrected  to  20°C  by  using  a  least  (squares 
regression  on  the  balance  housing  temperature  to  determine  the  temperature 
coefficient  for  the  measurements  (See  Figs.  11  and  12). 

This  incorporates  the  temperature  dependence  of  the  calibrations, 
which  otherwise  showed  no  significant  systematic  errors. 

The  mean  and  standard  errors(  o/v'rT  )  of  the  temperature  corrected 
measurements  are  given  in  Tables  3  and  4.  "N"  and  "A"  denote  null  and  active 
measurements  combined; 

Tables  5  and  6  give  the  %  drag  changes  observed.  The  errors  have  been 
calculated  from  the  standard  errors  combined  in  quadrature  and  expressed  as 
a  percentage  of  the  null  measurement:. 

It  can  be  seen  that  the  accuracy  of  these  measurements  is  much  better 
than  the  preliminary  data  of  ref.  1,  and  in  fact  at  high  flow  speeds,  the 
original  estimate  of  -  resolution  in  drag  change  has  been  bettered  to  ~  1/10% 

(ca.  5  mg  in  drag  force). 


6. 

* 

Fig.  13  shows  the  hull  measurements  in  the  form  D/P  (.«  C^)  vs  («  Re), 
showing  the  structure  of  the  flow  at  each  of  the  pitot  static  pressures. 

If  the  values  of  the  drag  coefficient  and  Reynolds  number  are  required 
explicitly,  comparison  can  be  mode  with  Fig.  10  .  The  data  from  the  first 
report  are  also  shown;  the  curve  is.  displaced  as  the  calibration  is  different. 

Another  series  of  measurements  at  very  low  flow  speeds,  using  the  deflection 
mode,  was  recently  initiated,  but  unfortunately  has  not  been  completed  due  to 
equipment  malfunctions.  It  is  hoped  to  finish  this  work  soon,  so  as  to  provide 
measurements  of  drag  changes  over  ns  wide  a  range  of  Reynolds  numbers  as 
possible,  and  to  improve  the  accuracy  of  the  1’  =  2  and  5  rnmllg  results. 


3.  TORSION  DISC  VISCOMETER 

3 . 1  Review  of  Work  by  Kestin  and  Shah 

Kestin  and  Shah^  made  an  experimental  investigation  of  changes  in 

the  apparent  viscosity  of  gases  produced  by  ionization.  They  used  a  torsion 

disc  viscometer  to  measure  the  apparent  viscosity,  and  radioactive  sources 

to  produce  ionization  of  the  gas  within  the  viscometer.  Two  approaches  were  used: 
137 

1.  a  Cs  Source  (  S-cmitter)  external  to  the  disc. 

210 

2.  a  Po  Source  (  a-emitter)  contained  within  the  disc. 

Measurements  were  made  on  the  following  gases  at  room  temperature  and  at  pressures 

-3 

ranging  from  1  to  1300  mmllg  (10  to  1,3  atm.):  Air  nitrogen  carbon  dioxide, 
oxygen,  helium*  neon,  argon,  krypton  and  xenon. 

Kestin  and  Shah  used  an  oscillating  (torsion)  disc  viscometer  principally 
because -of  the  simplicity  of  design,  and  also  because  the  results  obtained  are 
precise.  They  were  interested  to  determine  the  absolute  viscosities  as  well 
as  the  fractional  changes  induced  by  ionization,  and  so  they  had  to  make  elaborate 
calculations  to  obtain  the  former  from  the  measured  quantity  (the  logarithmic 
decrement,  or  damping,  of  torsional  oscillation).  The  dimensions  of  their  appara¬ 
tus  were  such  that  the  discs  were  of  3  to.  4  cm  diameter. 


The  results  obtained  by  Kestin  and  Shah  may  be  summarised  thus: 

1)  External  source  97.5  Curie  r  changes  in  apparent  viscosity  at 

1  atmosphere,  no  significant  change  (,<6i05%)  except  .  Kr  (+  0,22%}. h 
Xe  (+6.14%)  and  air  (-  0,25%) i 

2)  Internal  source  ca.  0.5  Curie  at  1  atmosphere,  mostly  detect 
significant  (  >0.5%)  changes in  the  range  -  1%  to  +  3%  at  1  mmHg, 
usually  increases  of  1%  to  2%  except  for  a  significant  descrease  of  7.6% 
in  the  case  of  air. 

3,2.  Introduction  to  the  Present  Torsion  Disc  Viscometer  Work 

It  was  decided  to  use  a  torsion  disc  viscometer,  again  for  reasons  of 
ease  of  construction  and  of  precision.  Since  we  did  not  require  to  know  the 
absolute  viscosity,  it  was  assumed  that  the  logarithmic  decrement  A  would 
be  nearly  enough  proportional  to  the  viscosity  over  the  likely,  range  of 
variation.  This  greatly  simplified  the  calculation.  The  theory  of  this 
form  of  viscometer  is  outlined  in  section  3,3  -  3.4. 

The  work  described  here  represents  principally  the  investigation  of 
a  suitable  design  for  a  viscometer,  with  a  number  of  preliminary  results  for 
atmospheric  air  at  room  temperature,  the  disc  being  irradiated  externally, 
by  (5  or  a  radiation. 

3.3.  Prototype  Design 

The  body  of  the  prototype  viscometer  was  made  from  dural.  The  disc 
itself  was  of  stainless  steel.  The  torsion  wire  was  held  by  two  pinchucks, 
one  fixed  in  the  torsion  head  and  one  in  a  brass  clamp  which  could  be  screwed 
down  onto  the  disc.  In  this  way  both  the  wire  arid  the  disc  could  be  changed. 
The  brass  clamp  was  fitted  with  a  small  plane  mirror  (aluminized  glass) .  The 
viscometer  was  fitted  with  two  'perspex'  windows,  to  allow  the  passage  of  a 
laser  beam  (He/Ne)  to  and  from  the  mirror;  this  beam  was  reflected  onto  a 
screen,  calibrated  in  millimetres,  about  two  metres  distant.  The  base  of  the 


•viscometer  ,  which  gave  access  to  the  disc,  was  a  push  fit  and,  like  .the  ' 
torsion  head,  sealed  with  a  rubber  'O’  ring.  (Wilson  seal).  The-i'hitial 
design  of  the  viscometer  is  ;shown  in  Figs..  14-16. 

Preliminary,  experiments  were  performed  using  various  torsion  wires  and 
a  disc  of-  diameter  63,5  mm  and  thickness  6.8  mm  (hom.)  Torsional  oscillations 
were  started  by  slowly  twisting  the  torsion  head  and  then  returning  it  to 
its  original  position,  as  nearly  as  possible.  The  procedure  was  found  to 
require  care  in  order  not  to  induce  too  much  vibration  in,  the  suspension; 
even  so,  it  was  usually  necessary  to  wait  about  10  minutes  for  the  grosser 
perturbations  to  die  away.  The  logarithmic  decrement  4  was  measured  by  con¬ 
verting  the  amplitudes  of  about  50  swings  of  the  tight-spot  on  the  screen  into 
the  angular  amplitude  of  oscillations  of  the  disc.  A  graph  (Fig  17)  shows 
typical  results.  The  results  of  the  preliminary  experiments  are  given  in 

**v* 

Table  7. 

At  this  stage  some  modifications  were  made.  It  was  thought  advisable 
to  have  a  means  of  levelling  the  viscometer  accurately.  The  disc  itself 
was  used  as  a  plumbr.line.  A  'perspex'  insert  2  cm  long  was  put  into  the 
shaft  of  the  viscometer.  The  external  section  of  this  insert  was  square,  with 
a  thin  vertical  line  scratched  and  inked,  centrally  on- each  face.  The 
viscometer  was  stood  on  a  levelling  plate,  which  was  adjusted  until  the 
torsion  wire  was  coplanar  with  each  of  the  two  pairs  of  opposite  lines. 

The  viscometer  was  next  fitted  with  a  gas  port  -  a  short  length  of  pipe 
sealed  and  welded  to  the  upper  surface  of  the  body.  It  was  connected  to  a- 
manometer  and  partially  evacuated:  in  this  way  it  was  possible  to  locate  and 
seal  any  leaks  with  "Araldite", 

In  order  to  reduce  the  labour  involved  in  calcula  i ug  A  from  the  data 

a  computer  program  was  written,  and  several  experiments  were  performed  to 

find  how  accurately  and  how  reproducibly  A  could  be  measured.  Preliminary 

attempts  were  made  to  find  the  effect  on  A  of  placing  radioactive  plates, 

.  .  147 

containing  Pm  (a  (3  emitter)  at  an  activity  of  about  150  mCi  on  the  base 


of  the  viscometer,  and  -then  comparing  the  results  with  those  from  "dummy"  ' 
(inactive)  plates.  This  involved  removing  the  base-  from  the  viscometer 
between  experiments,  which  was  unsatisfactory,  since  it  disturbed  the 
torsion  wire.  Experiments  were  also  tried  with  the  base  removed  and  Che 
viscometer  positioned  above  the  radioactive  plates.  The  results  obtained 
are  contained  in  table.  8, 

It  was  found  that  the  values  of  A  recorded  in  these  experiments  varied 
by  about  6%.  The  reason  for  this  was  not  known,  but  it  was  conjectured 
that  it  was  due  to  vibration.  In  an  attempt  to  reduce  the  effect  a  heavier 
disc  (3.2  mm  (nom.)  thickness)  was  used,  and  a  stand  was  designed  to  isolate 
the  viscometer  from  vibration  and  to  aid  temperature  control.  Finally  the 
base  of  the  stand  was  made  a  screw  fit  so  that  it  could  be  removed  and  replaced 
with  the  minimum  of  disturbance  (Fig.  18). 

Further  experiments  were  performed  with  the  heavier  disc,  using  38  gauge 

(0.15  nun)  Cu/Be  wire.  The  reproducibility  was  a  little  better  than  before. 

58 

An  encapsulated  source  of  Co  (y  emitter;  about  0.3  mCi)  was  compared  with 
a  dummy  cut  from  dural.  A  thermocouple  was  led  through  the  gas  port  in 
order  to  assess  any  correlation  between  A  and  temperature.  Also,  experiments 
were  performed  on  the  computer  to  asses  the  effect  on  the  measured  value  of  A 
of  misreading  the  zero  position  on  the  scale  or  of  incorrectly  measuring  the 
distance  from  the  viscometer  to  the  scale.  These  experiments  were  made  by  using 
one  set  of  data,  and  varying  these  parameters  each  time  they  were  put  through 
the  computer.  The  results  of  these  experiments  are  contained  in  Table  9  and  10. 
A  significant  correlation  was  found  between  A  and  temperature  and  this  is 
shown  in  Fig.  19  and  corrected  for  in  the  final  results. 


V-  ,"-v l  .^v  >  A!  v 
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3.4  Theory  of  :the  Torsion-Disc.  Viscometer 


A  , disc  of  radius  a  experiences  a  couple  -  p0' 
when  it  is  '.twisted:  through  an  angle  9  in  :its  own. 
plane,  about  its  centre.  It  lies  between  two 
parallel,  fixed  discs  each  at  a  distance,  d.  Suppose 
this  disc  has  an  angular  velocity  w  and  consider 
the  force  acting  on  the  element  of  one  surface  lying 
between  circles  of  radii  r  ,and  r  +  dr. 


Velocity  at  surface  of  disc  =»  s>r 

velocity  at  fixed  disc  =  0  (no-slip  conditions) 

, * .  velocity  gradient  « 

dF  =  q-^~  *  2jrrdr 

'**■-♦*  a 

Thus  the  torque  acting  on  the  element  is 

dG  =  2H-r3dr 

,  and  the  total  torque  acting  on  both  sides  of  the  disc  is 
Thus  we  can  write  the  equation  of  motion  of  the  disc: 


10  + 


0  +  p0  =  0 


where  I  is  the  moment  of  inertia.  The  solution  is  a  dampled  simple 


harmonic  motion 


where 


0  =  0Qe  sin(wt  +<{>)  (for  light  damping) 


a 


A  =  2iia /to 


2  4 

TT  a 

— — -  dp)  2.n. 

a 


In  other  words  A  is  proportional  to  h, In  practice,  the  damping  term  is  of 


the  form  ^ — j-+k|0  >  anc^  so  =  cr>  +  4^,  but  the  term.  Aq  (arising  from 
mechanical  defects  in  the  wire,  and  so  oh)  is  sually  very  small  '(typically 


5  x  10~5). 


More  importantly,  "edge  effects"  should  be  taken  into  account,  which 


depend  upon  the  finite  thickness  of  the  disc.  However.,  these  merely  have 


the  effect  of  altering  the  calibration  of  the  instrument  with  the  viscosity 


and  density  of  the  gas  contained,  and  since  these  change  by  .the  order  of  1%, 


■  this  is  an  error  in  the  measured  change  of  viscosity  which  will  be  due  to 


edge  effects.  On  the  other  hand,  if  an  absolute  measure  of  viscosity  is  requirl 


it  will  be  necessary  to  produce  empirical  claibration  curves. 


4.  SUMMARY  OF  EXPERIMENTAL  RESULTS 


p 

The  results  of  the  drag  balance  measurements  are  summarised  in  Figs.  20  anjj 


21.  It  can  be  seen  that  for  both  sets  of  measurements  the  drag  changes  for 


P  ^lOmmllg  (ca.  50  m/s)  are  entirely  consistent  with  no  change  in  drag  force 


on  irradiation,  and  that  any  change  for  the  higher  flow  velocities  is  less  than] 


ca.  0.2%.  For  P  =  2  mmHg  there  is  an  apparent  drag  reduction  of  a  few  %, 


most  significant  for  the  turbulent  flow.  This  measurement  is  still,  however, 


less  than  2  standard  deviations  from  no  change.  The  increase  for  P  »  5  mmHg 


(laminar  flow)  in  conjunction  with  the  P  =  2  mmHg  decrease  could  be  taken  as 


evidence  for  a  change  in  the  structure  of  the  developing  turbulent  flow 


(see  Fig.  13). 


The  torsion  disc  measurement  with  the  'Pm  plates  gave  a  change  of  vis-  ): 


+  ‘  58  +  ‘ 

cosity  of  +  0.08  -  0.26%  while  for  Co  the  change  was  0.0  -  1.2%,  both  consist* 


with  no  change. 


DISCUSSION  AND  CONCLUSION 

As  this  report  describes,  we  have  now  completed  the  construction  and 
evaluation  of  a  drag  balance  and  a  torsion...  disc  viscometer*  Both  are 
instruments  of  high  precision  which  allow  us  to  study  changes  of  drag  ;and 
viscosity  in  gases  to  high  accuracies.  To  date  they;  have  been  used  to 
investigate  aiiy  changes  caused  by  radioactive  irradiation  at  atmospheric 
pressure. 

(A) 

Both  the  theoretical  and  experimental  work  of  Kestin  and  Shah  suggest 
that  only  a  very,  small  drag  change  is  to  be  expected  at  atmospheric  pressure 
with  the  radiation  intensities  used.  This  is  in  agreement  with  calculations 
using  the  Chapman  Enskog  theory  (see  Ref.  5)  which  suggests  that  the  size  of  the 
drag  (or  viscosity)  change  is  of 'the  order  of  the  ionisation  intensity.  The 
present  experimental  work  seems  to  be  in  agreement  with  these  results,  although 
there  is  a  possibility  of  an  anomalous  drag  reduction  at  low  flow  speeds. 

We  now  plan  to  extend  our  research  to  lower  speed  measurements  and  also  to 
begin  work  at  pressures  less  than  atmospheric.  These  experiments  should  prove 
exciting  since  there  is  evidence  of  a  viscosity  decrease  for  air,  caused  by 
radioactive  bombardment  at  lower  pressures  (ref.  4).  The  research  will  be 
directed  at  finding  the  optimum  conditions  for  any  decrease  of  drag. 

Experiments  are  also  planned  to  establish  the  basic  mechanisms  involved 
in  any  drag  reduction.  The  three  possibilities  which  will  bo  examined  initially 
•are  that  it  is  due  to  (i)  ionization  (ii)  thermal  input  at  a  critical  part  of 
the  boundary  layer  (iii)  molecular  clustering.  The  apparatus  described  in  this 
report  has  the  flexibility  to  test  these  various  possible  mechanisms  and  with 
techniques  already  available  in  the  laboratory,  such  as  high-speed  photography, 
flow  visualization  and  mass  spectrometry,  it  should  be  possible  to  make 
a  full  investigation  of  drag  reduction  processes. 
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TABLE  $ 


%  DRAG 

CHANGES  ON  IRRADIATION  (FULLY  TURBULENT' FI OW ) 

(maximum  accuracy  quoted  nearest  0.05) 

P/mmHg 

A1 

A2 

A. 

2' 

“1 #4  +  1,1 

-1.6  +  1.1 

-1.5  +  0.8 

5 

*0. 75+  0.5 

+0.2  +  0.35 

-0.3  +  0.-35 

10 

+0,35^+  O.45 

+0.5  +  0.35 

+0,4  +  0,-3 

'  15 

.  +O.05+  0,2 

-0.1  +  0,15 

-0.05+  0.15 

50 

-0.25+  0.2 

-0.15+  0.15 

-0.2  +  0.15 

50 

+0,1  +  0.2 

6.0  +  0,3 

+0.05+0.2 

TABLE  6 

• 

%  DRAG  ( 

CHANGES  ON  IRRADIATION 

(LAMINAR/TRANSITI6NAI.  FLOW! 

(maximum  accuracy  quoted  nearest  0,05) 

"  *  ' 

P/mmflg 

'.Al 

A  2 

A 

2 

-0.4  +  2.8 

-5.0  +  1.9 

-1 .8  +  1,9 

5 

+1.7  +.  1.7 

+2.0  +  1,2 

+1*9  +  1.1 

10 

-0.6  +  1.1 

+0.9  +  0.8  ' 

+0.15+  0.8 

15 

-0.6  +  0,65 

-0.2  +  0.75 

-0.4  +  0.55 

50 

+0.05+  0.3 

-O.O5+  0.25 

0.0  +  0.2 

50 

+0.1  +  o.l 

-0.2  +  0.1 

-0.05+  0.1 

TABLE,  7 


Preliminary  experiments 


wire  .  dia/rom  .  ,  ...  ...  A, 


steel (piano) 

0.127 

0.0088  , 

0.0058 

98Cu/2Be 

0.051 

0.030 

it 

0.149 

6.048  + 

0.0006 

'  table. 8  ■ 

.  .  147 

Experiments  with  Pm 


Axio6 

Base  prized  off 

between  experiments 

dummy 

4999  +  9 

active 

5053  +‘11 

6 A  =0.995  +  0.003 

dummy 

5156  +  9 

A 

Base  removed; viscometer  propped 

dummy 

4693  +  17 

active 

4824  +  35 

6A  =1.030  +  0.008 

dummy 

4676  +  17 

A 

dummy 

5078  +  63 

active 

5130  +  31' 

6 A  =0.995  +  0.010 

dummy 

5229  +  49 

A 

Base  prized  off 

between  experiments 

;heavy  disc 

dummy 

2866  +  23 

active 

2920  +  18 

SA  =1.016  +  0.008 

dummy 

2884  +  22 

A 

Combined  result:  Change  in  viscosity  on  irradiation  from 

+(0.08  +  0.26)% 


•TABLE.  9. 


Sensitivity  to  misreading 

z  =  centre  of  oscillation  on  scale  ,  cm  -. 

1  =  distance  from  viscometer  to  scale  ,  cm  . 


1 

AxlO6 

>9.0 

2-)l 

2500 

59.1 

257 

2616 

39.9 

257 

2716 

40.0 

257 

2742 

41.0 

257 

5025 

40.0 

252 

2742 

40.0 

242 

2742 

1  accurate  to  £  5  cm  =>  A  accurate  to  <  0.05% 
z  accurate  to  _+  1  mm  =>  A  accurate  to~l  % 
z  accurate  to  +  1  cm  A  accurate  to —9  % 


TABLE  10 


Expt .code 

AxlO6 

T  C 

D00  ' 

3 030 

"11.7 

dummy 

Dll 

2991 

10.6 

dummy 

D12 

2931  . 

12.0 

dummy 

D13 

2986 

15.0 

dummy 

D14 

2947 

16.7 

dummy 

Di5 

2857 

19.2 

dummy 

Di6 

2956 

20.7 

dummy 

D21 

2665 

20.2 

active 

D22 

2877 

20.5 

active 

D23 

2716 

21.0 

active 

D24 

2772 

21.6 

active 

D25 

2934 

22.0 

active 

D26 

2904 

22.0 

active 

D27 

2656 

21.9 

active 

D31 

2831 

17.9 

dummy  ■ 

D32 

2742 

l8.6 

dummy 

D33 

2812 

18.8 

dummy 

D34 

2945 

20.6 

dummy 

D41 

2959 

l8,6 

active 

D42 

3025 

18.6 

active 

D43 

2796 

16.8 

active 

D44 

2920 

16.4 

active 

D6l 

2897 

15.3 

dummy 

D62 

2907 

■17.1 

dummy 

D63 

2826 

21.9 

dummy 

D64 

v  2837 

21.3 

dummy 

D71 

2932 

18.9 

active 

D?2 

2897 

20.0 

active 

D73 

2796 

21.4 

active 

D74 

3023 

22.1 

active 

A  =  2900 

+  21 

dummy 

A  =  2858 

+  31 

active 

uncorrecfced  for  temperature 


Regression  of  A  on  T  : 

A  =  3141  -  14.0  T 
r2=  0.37 

values  of  A  corrected  to  20  C  : 
&  =  2860  +  16  dummy 

A =  2860  *  30  active 


Change  of  viscosity  on  irradiation  from 

(0.0  +  1.2)% 


